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Abstract 

The bilinear combination of Dirac spinors u{pi,ni)u{p2,n2) is expressed in terms of Lorentz vectors 
in an explicit covariant form. The fact that the obtained expression involves only one auxiliary vector 
makes it very convenient for analytical computations with REDUCE (or FORM) package in the helicity 
formalism. The other advantage of the proposed formulas is that they apply to massive fermions as 
well as to massless fermions. The proposed approach is employed for the computation of one-loop 
Feynman diagrams and it is demonstrated that it considerably reduces the time of computations. 


1 Introduction 

Over many years, expressions of the type 

w{pi, ni)w{p2, n2) = S + (1) 

have received considerable study [l]-[6]. They have been extensively used both in calcula¬ 
tions of multiparticle helicity amplitudes and in the analysis of polarization phenomena. The 
method of computation of multiparticle helicity amplitudes (hereafter referred to as the spinor 
formalism) is based on the representation of the polarization vectors of gauge bosons in terms of 
spinors. It has been successfully applied to the computations of tree amplitudes of the reactions 
involving massless particles [2] and has also been extended to the massive case [7, 8]. However, 
a computation of a loop diagram with this method presents a challenge. 

Another approach is to express the relevant products of spinors in terms of Lorentz vectors 
(the method of calculations with the use of such expressions is referred to as the helicity 
formalism). However, the expressions for S, etc. obtained in [3, 4] are very cumbersome and 
involve singularities that hinder employing these expressions in computations. One can use a 
particular basis to avoid the singularities; in this case, there emerges an additional auxiliary 
vector [5, 6]. 

In this work, we discuss properties of these auxiliary vectors and demonstrate that an ex¬ 
pression of these vectors in terms of momenta simplihes the computations and makes it possible 
to employ the helicity formalism for an analitycal computation of the helicity amplitudes in the 
one-loop approximation. 

^E-mail: rogalyov@mx.ihep.su 


1 



2 Basic Formulas 


In order to compute a complete set of the helicity amplitudes of some process involving 2 
fermions, it is sufficient to consider the case when the momenta pi and p2 and the spin vectors 
±A^i and ±iV2 of these fermoins lie in the one 2 -plane; in this case, polarization vectors of the 
fermions can be expressed in terms of the momenta by the formulas 


7Vi = 


1 

A 


f Pl • P2 

\ mi 


Pi - mip2 , 


N 2 


1 

A 


(^-m2Pi 


+ 


Pi ■P2 
m2 



( 2 ) 


where m^ = pf, m\ = p^, and A = y (pi ■ ^2)^ — mlml- This assertion stems from the fact 
that an arbitrary polarization state of a fermion is a quantum-mechanical superposition of the 
states \N) and | — N), where \N) is the state with the spin vector directed along some vector 
N and | — N) is the state with the spin vector directed oppositely to N. For this reason, we 
restrict our attention to the basis composed of the fermion states |pi,iVi), |pi, —Ni), |p2,-^2), 
and \p2, —N2) (or, in another notation, u{pi, Ni), u{pi, —Ni), u{p2, N2), u{p2, —N2)), where the 
vectors Ni and N2 are defined by ( 2 ); this basis is named ‘diagonal spin basis’ [ 5 , 6 ]. 



Figure 1 : Momenta Pi,P2 smd vectors of spin Ni, N2 lie in the one 2 -plane 


In this basis, the combination of Dirac spinors u{pi, Ni)u{p2, N2), where pi and p2 are the 
momenta and Ni and N2 are the polarization vectors, has the form [ 9 ] 


/l±7® It 7^ '' 

uu{±, ±) = ji —--j2—-h mik2 - ^2^1 


V 2 ’ 


( 3 ) 
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mm(±, t) = {jiki + mi)k2 —^-h {j2k2 + m2)ki —^— 

where ««(+, —) = m(pi, ni=+iVi)-u(p2,''^2=—-^2) etc., 


( 4 ) 


Ji = 


VPi ■ P2 + mim2 + VPi ■ P2 - ■nrim2 

71 ^ 


VPi ■ P2 + rnim2 - VPi ' P 2 - mima 
J2 =-^-, ( 5 ) 
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and the light-like vectors /ci, ^2, and ui± are given by 


2A 


mi , 


m2 


ki = 7 ^[ jiPi -J2P2 , k2 = - j2Pi + jiP2 , ki-k2 = 


2A 


m2 . 


mi 


( 6 ) 


= 


7^3 


i± 77 r a 

7gK2— h k2qki —-— 


( 7 ) 


where 

7^3 = ■ ^2(2 ki-qk 2 -q-q-qki-k 2 ) (8) 

and q is an arbitrary vector such that A3 7^ 0 ^. To simplify computations, one should choose 
vector q such that either g ■ g = 0 or g ■ /ci = g ■ A;2 = 0 . In spite of the presence of an arbitrary 
vector g in ( 7 ), a;+ and are almost independent of g: they only acquire a phase as the vector 
g varies. 

The explicit expressions ( 7 ) for uj± allow to compute the complete set of the helicity am¬ 
plitudes using only momenta of the particles, that is, without resort to specihc polarization 
vectors [ 2 ] or ‘fundamental’ spinors [ 7 ]. Such computation is possible due to the arbitrariness 
in a choice of the ‘gauge’ vector g. To explain this in more detail, let us consider a diagram 
involving a fermion line of the type T = u{p2,n2)Ou{pi,ni), where O is something like T in 
formula (10) below. Since T = Tr Ou{pi,ni)u{p2in2)i it can be evaluated with the formulas 
( 3 )-( 7 ); this being so, the momentum of any particle other than a fermion of momentum pi or 
P2 can be inserted in ( 7 ) for vector g. 

The only exception from this rule is the decay fi /27, where fi and /2 are heavy and light 
fermions and 7 is the photon (or gauge boson). In this case, there is no momentum independent of the 
fermion momenta and, therefore, vector q with the required properties does not exist. In the analysis 
of such reactions, the vectors uj± should be identified with the polarization vectors of the photon. 


3 An Example of Computations 

Our approach is based on the representation of spinors in terms of vectors, its application 
to loop computations is straightforward. Therein lies the difference between the proposed 
approach and that by Ballestrero et al. [ 7 ], in which all relevant vectors should be expressed in 
terms of spinors. We illustrate the proposed method by computing the imaginary part of the 
one-loop diagram for the decay K{pk) —^ p,{k)h'{k')'y{q) (see Fig. 2 ). The kinematical variables 
here are defined as follows: 

2 pK ■ g 2 pK-k mj 

A = A^; C=( 1 -A)( 1 -i)-p, 

_ X 

^Here and below, 7® = i7°7^7^7^ and case, Tr7^7^7“7^7^ = and the operator 

g{l,a) in the RELUCEpackage (symbol I marks the fermion line) should be identified with —7®. 
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where is the lepton mass and is the kaon mass. 

We employ the Cutkosky rules [ 10 ] to replace cut propagators by the 6 functions. Thus we 
obtain the expression for the imaginary part of the amplitude, 


M" = ^G,eV:, n(k')(l+',^)T(r,k,k',q,e)v(k), ( 10 ) 

2 n J Z[r-q,r-k) 


where T = 8(fc + g)(l — — q — k)'y’^{m£ — q — f)e{m£ — f)'y^, Z is the product of the 

remaining (uncut) propagators, r is the loop momentum, spinors u and v describe the neutrino 
and muon, and e = e^(±)7^. Here 

^ 2M^a/AC ^ ~ - (1 - P - T • (11) 

is the polarization vector of the photon. The factor ^GpeV*^ in formula ( 10 ) should be thought 
of as an effective coupling constant. 



Figure 2 : Diagram giving a contribution to the amplitude of the decay K —/iz/7. 


Our method is based on computation of the helicity amplitudes.^ Making use of the formulas 


v^{k, -N)u^{k‘ 
v^{k,N)uu{k')- 


{k — m£)k' 


2 MKy/l — X — p 
M|;(l — X — p) — mek' 


2Mk\^1 — X — p 


a+1% 

i"-(l + 7 ^), 


( 12 ) 


^Here there are four helicity amplitudes: the muon helicity (we select the reference frame comoving with the 
center of mass of the muon and neutrino) takes one of the two values ±1/2; photon helicity in any reference 
frame is ±1. 
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where 


Cc;_ = 


2M]^pxX 

i ^ 

1 — X — p 


(13) 


and 


V 2 

2Mlx^ 


{hqk'il + 7 ^) + k'q^k(l 





{l-x- p)k^ - 2pkl 

mi{l — X — p) 


(14) 


we express the quantity u{k')T{r,k,k',q,e)v{k) (which is nothing but the respective helicity 
amplitude) in the integrand in (10) in terms of the scalar products of k, k', q, and r for each 
hxed polarization state of the muon and photon and only then perform integration with 
respect to r. The two integrations are trivial due to the <5 functions, the integtrand takes the 
form of a rational function F = F{x, X, p, v, cos0): 


dr6{r‘^ — mf)5 (^{k + q — r)^) F{r) 


x{l — X) 
8 r 


dv 


r27T 


/-I 


d(f)F^ 


(15) 


where 

= ill _ 2p _ 

x\l-XY Ml x{l-X) ’ ^ ^ 

and the azimuth angle 0 specihes the direction of the projection of the vector r onto the plane 
orthogonal to the 4-vectors k and q. The function F{x, X, p, Mk/Mt^,v, coscf)) is a polynomial 
in COS0, therefore, integration with respect to 0 is trivial; this being so, the integrand in the 
integral with respect to v has the form A{y) + Bj{Cv F D), where A{v) is a polynomial in v 
and B, d, D are independent of v. The computation of the diagram in Fig. 2 is made with the 
REDUCE package [11]. This diagram is calculated exactly, no approximation is used [12]. A 
computation in the helicity formalism is by an order of magnitude faster than a computation 
with the traditional method. 

Now we present the contribution of the diagram in Fig. 2 to the imaginary part of the 
helicity amplitudes as it was calculated using the above procedure: 


The variables used below are defined as follows: 
let w=l-lmbd; 
let b~2=x*lmbd; 
let c~2=x*(l-lmbd); 
let az"2=l-rho-x; 
let d~2=c~2*az''2-rho*b''2; 


7o7o7o Results of computations: 

GammaMinusMuMinus := (16*rho**3*me*mk**4*w*x**2*( - w + 1) 

+ 16*rho**2*me*mk**4*w*x**2*( - w**2*x + w*x + w - 1) 
+ 4*rho*me*mk**4*w**2*x**4*(w**2 - 1) 

+ 8*me*mk**4*w**3*x**4*(w*x - w - x + 1)) 

/(rho*az*b*d + az*b*d*w*x)$ 


GammaMinusMuPlus := (16*rho**2*mk**5*w*x**2*(w - 1) 
+ 4*rho*mk**5*w**2*x**3*(w - 1)) 

/(rho*az + az*w*x)$ 
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GammaPlusMuMinus := (16*rho**3*me*mk**4*w*x**2*( - w + 1) 

+ 8*rho**2*me*mk**4*w*x**2*( - 2*w**2*x - w*x + 2*w + 3*x - 2) 
+ 4*rho*me*mk**4*w**2*x**3*(w**2*x - 9*w*x + 6*w + 8*x - 6) 

+ 4*me*mk**4*w**3*x**4*( - w*x + w + x - 1)) 

/(rho*az*b*d + az*b*d*w*x) 

+ ln(w*x/rho+l) * (8*rho*b*d*me*nik**4)/az$ 


GammaPlusMuPlus := (8*rho**2*mk**5*w*x**2*(w + 1) 

+ 4*rho*mk**5*w**2*x**2*( - w*x + 6*x - 4) 

+ 4*mk**5*w**3*x**3*(x - 1)) 

/(rho*az + az*w*x) 

+ ln(w*x/rho+l) * ( - 8*rho**2*nik**5*x + 8*rho*mk**5*w*x*( - x + l))/az$, 


where Imbd = 1 — A, rho = p, me = mk = the variable GammaPlusMuMinus is the 
imaginary part of the amplitude for the photon helicity +1 amd the muon helicity —1/2 (in 
the reference frame comoving with the center of mass of the lepton pair) etc. Here we avoid 
cumbersome denominators by expressing all vectors in terms of the light-like vectors q, 


ki = 


1 

Mk X - p 



1 - 


P _ 

X — p 



and k 2 


k' 

Mk X - p 


(17) 


from the outset. 


Note added 

The imaginary part of the above diagram is needed for the determination of the transverse 
component of the muon spin, which may provide an indicator of CP violation. It should 
be noted that, using the relations similar to (3) and (4), the average value of the transverse 
component of the muon spin*^ in the reaction under consideration can be expressed in terms of 
the helicity amplitudes as follows: 

^ ^ (^M'_ _M'L + - M'_ +M" _ +Af'+ _Ml + - +M" _) , ( 18 ) 

where Af is the normalization factor, = J2ij=± and and AA'^^ are the real and 

imaginary parts of the helicity amplitudes, AAr^s = JA'rs + ^■^'1 s (indices (r, s = ±) mark the 
photon and muon helicities). 

Appendix. Properties of the Vectors 

Any ordered pair of real light-like vectors ki and k 2 '. ki ■ k 2 = \ determines the other pair 
of the complex zero-norm vectors uj+ and a;_, which are complex-conjugated to each other and 
satisfy the relations 


■ ki = ■ k2 = 0, ■ ki = ■ k2 = 0, u;_ ■ a;+ = —1. (19) 

^It measures a half of the muon polarization 
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Now, let the tensors 

kioj’^ — oj^ki and kt^oo'^ — oo^k^ 

be self-dual: 

kiu+{l-'y^) = k2(2}-{l-'y^) = 0, ( 20 ) 

whereas the tensors 

^2^^+ — ^+k 2 and kiu'^ — u^ki 

be anti-self-dual: 

fc2^h+(l + 7^) ~ ki(jj-{\ 7^) = 0. (21) 

The vectors ujj^ and uj^ are determined by the pair ki, k 2 and the requirements (19) 
and (20)—(21) up to a phase. 

• The vectors {oo-) can be interpreted as the polarization vectors of the photon of 
momentum ki and positive (negative) helicity (in the gauge A;2 ■ ^ = 0) 

The vectors uj± can be represented in the explicit form® [6, 9] 

= -^{k,-k2q^-q-k2k>(-q-k,k!)±ie^^^^^^), ( 22 ) 


where A3 is dehned in (8). 

It should also be mentioned that, in the framework of the Cartan theory of spinors [13], a 
Dirac spinor in 4-dimensional space can be interpreted as an isotropic 2-plane in 5-dimensional 
pseudo-Euclidean space (4 dimensions in it are associated with the Minkowski space-time and 
the 5th dimension is associated with 7® giving the signature (-I-, -|-, —, —,—)). This being so, 
such plane for the spinor® u(p,n) is nothing but the linear span of the zero-norm 5-vectors (1, n) 
and (0,a;+), where the 4-vector is associated with the vectors 

h = ^ and k 2 = ]r (— -n^ 

2 \m J 2 \m ) 

according to the above procedure. 
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